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Abstract 

We determine viscosity corrections to tlie retarded, advanced and symmetric gluon self energies 
and to the static propagator in the weak-coupling "hard loop" approximation to high-temperature 
QCD. We apply these results to calculate the imaginary part of the heavy-quark potential which is 
found to be smaller (in magnitude) than at vanishing viscosity. This implies a smaller decay width 
of quarkonium bound states in an anisotropic plasma. 
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Introduction: The retarded, advanced and symmetric gluon self energy in an equilibrated, 
weakly-coupled high-temperature plasma can be calculated within the "hard thermal loop" 
effective theory |2!]. However, the standard expressions receive corrections if the plasma 
undergoes (anisotropic) expansion and if its shear viscosity is non-zero. In this letter, we 
provide explicit expressions for the leading viscous corrections to the retarded, advanced 
and symmetric gauge-boson self energies and to the corresponding (static) resummed prop- 
agators. As an application of particular interest, we finally determine the corrections to 
the imaginary part of the heavy-quark (static) potential in a QCD plasma. This imaginary 
part provides a (contribution to the) width F of quarkonium bound states [l, 0, [sj which 
in turn determines their dissociation temperature: dissociation is expected to occur when 
the binding energy decreases, with increasing temperature, to ~ F j^, 0]. It is therefore 
interesting to determine the effects of non-zero shear viscosity on the imaginary part of the 
potential. 

Corrections to the thermal distributions in an anisotropic plasma: We consider a hot 
QCD plasma which, due to expansion and non-zero viscosity, exhibits a local anisotropy in 
momentum space. The phase-space distribution of particles is given by 



/(P) = /isoU/p'+^(P-n)' (1) 



/iso(p) 



l-^^^^(l±/i.o(p)) 



(2) 



Thus, /(p) is obtained from an isotropic distribution /iso(|p|) by removing particles with 
a large momentum component along n, the direction of anisotropy [8]. We shall restrict 
ourselves here to a plasma close to equilibrium and so /iso(p) is a thermal ideal-gas distribu- 
tion equal to either a Bose distribution nsij)) or to a Fermi distribution urip), respectively. 
Equation ([21) follows from an expansion in the anisotropy parameter ^. The correction Sf 
to the equilibrium distribution exhibits precisely the structure expected from viscous hydro- 
dynamics for a fluid element expanding one-dimensionally along the direction n, provided 
we identify jof 

5^^^. (3) 

Its 

Note that this expression only holds true in the Navier-Stokes limit. In the general case. 



one can relate ^ to the shear tensor [10|. Here, 1/r denotes the expansion rate of the 
fluid element and 77/5 is the ratio of shear viscosity to entropy density. As usual in viscous 
hydrodynamics, the temperature T as well as the entropy density s appearing in Eqs. ([2]) 
and ([3]) are defined at equilibrium; all viscous corrections are accounted for explicitly by 6f 
in Eq. ^. 

Propagators in the Keldysh real-time formalism: We shall calculate the finite-temperature 



self energies and propagators using the well-known Keldysh real time formalism The 
propagators are then 2x2 matrices such as 

OiP)-(^'J^)- HP-' - ™=) ( 3^ - ) (4) 



for a scalar field and 



^ ^ V -0(Po) + fp fp 



(5) 



for a Dirac field. We use the following notation: P = {po,p), p = |p|. In equilibrium, 
the distribution functions fs and fp correspond to Bose or Fermi distribution functions, 
respectively. Away from equilibrium they need to be replaced by the corresponding non- 
equilibrium distributions from viscous hydrodynamics. It should be noted that Eqs. (jlj) and 
dS]) are "bare" propagators, the hard loop resummation has yet to be performed. 

The retarded, advanced and symmetric propagators can be obtained from the Keldysh 
representation (which satisfies Du — D12 — D21 + D22 = 0) via 

Dn = Dn - D,2 , Da = Dn - D21 , Dp = Dn + D22 . (6) 

In momentum space, the explicit expressions are 

Dr{P) ^ 



P"^ — m? + i sgn(po)e ' 
Da{P) = 2 '■ 7 — T' 

Df{P) = -2m (1 + 2/b) 5(P' - m^) (7) 



for scalar bosons and 

Sr{P) 
Sa{P) 



f + m 



P"^ ~ m? + i sgn(po)e 
f + m 



P'^ — m? — i sgn(po)e ' 
Sp{P) = -27Ti if + m) {1 - 2fF) S{K^ - m^) (8) 

for fermions, respectively. 

In the real time formalism, similar relations hold for the self energies: 

Hn + + Hai + Has = (9) 

and 

= Uu + , = nn + 1I21 , = nn + 1I22 . (10) 

Resummed photon ( gluon ) propagator: The resummed photon propagator can be deter- 
mined from the Dyson-Schwinger equation 

iD* = iD + iD {-m) iD\ (11) 

where the propagators and self energy are 2x2 matrices. D* indicates a resummed propa- 
gator and D is the bare propagator. 

Using the identities ([6]) for the bare and resummed propagators and ffTOl) for the self 
energies, it is easy to show that 

D*R = Dr + DrUrD*r. (12) 
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A similar expression holds for the advanced propagator. The resummed symmetric propa- 
gator satisfies 

D*p = Df + DRYiRD*F + DfHaD^a + ^hH^D*^ . (13) 

Using Df{P) = (1 + 2/b) sgn(po) [Dr{P) — Da{P)] (this equation is true even in the non- 
equilibrium case lH), the solution for D*f can be expressed in the form 



D*FiP) = (1 + 2/b) sgn(po) [D*RiP) - D* a{P)] 

+D*r{P) {Uf{P) - (1 + sgn(po) [IIr{P) - TIa{P)]} D*a{P) . (14) 

In equilibrium we have Hf^P) = [1 + 2nB(po)] sgn(po) \^r{,P) — ^a{P)] and as a result the 
second term in Eq. f|T^ vanishes. However, this is no longer true out of equilibrium. 

For the static potential we only require the temporal component of the gluon propagator. 
At leading order in ^ (or in the shear viscosity) the calculation of D*^^ is easiest in Coulomb 
gauge. The temporal component decouples from the other components which simplifies 
the calculation significantly. In this gauge (with the gauge parameter r/ = 0) the bare or 
resummed propagators satisfy 

K'-D^' = ,i = 1,2,3 (15) 

as a consequence of diA^ = 0. In the isotropic case this reduces to = 0. We can then 
write Eq. fll2l) as 

D%) = D^^ + D^n^(o)I)*^(o) . (16) 

Here, L denotes the temporal component, = D^^. A similar relation holds for the 
advanced propagator. Eq. f|T6l) no longer holds for the anisotropic system due to breaking of 
isotropy in momentum space. However, for small anisotropy, we can expand the resummed 
propagators and self energies in ^: 

D* = z}(*o) + ^D^^ + 0{e), n = H(o) + en(i) + 0{e) . (17) 

The propagators to order V^q^ (either retarded, advanced or symmetric) satisfy the 
equilibrium relations mentioned above. For the linear term of order ^, 

D*R{i) = -Di? H^(i)D*i{(o) + -Dij Hi{(o)-D*ij(i) . (18) 

For the temporal component, we have 

D%) = {Dr)'^ (n«a)),.(D*K(o))^° + {DRf- (H^(o)),.(D*K(i))^° . (19) 

Since {Dr)^' = 0, {D*r^o))^' = and (Hr(o))°^ ~ P\ according to Eq. it follows that^ 

-^*R(l) = ^R^R{l)^*RiO) + ^R^RiO)^*Ril) ■ (^O) 



^ The temporal component of the retarded propagator to order -0*^(2), fails to satisfy such a relation. 
In fact, it includes a product of {Dr)°'^ (nj^(i))^^(D*fl(i))''°, but (nfl(i))°* is not proportional to P\ As 
a result, n^(]^)I?*/{(i) doesn't give zero automatically. This term will therefore depend on the spatial 
components of the self-energy and propagator which makes the calculation more complicated. 



4 



Again, a similar relation holds for the advanced propagator. For the symmetric propagator, 
finally. 



D 



F(l) 



;i + 2/^(0)) sgn(po) [D%,^{P) - D%^{P)] 



+ 2fBa)Sgn{po)[D%){P)-D* 

+ D%^{P) {n^(i)(P) - (1 + 2/^(0)) sgn(po) [n^(i)(P) - ni(,)(P)] 

- 2/5(,)Sgn(po)[n^(o)(i^)-ni(o)(i')]}^*A(o)(^) • 



(21) 



We now proceed to calculate explicitly the photon/gluon self energies from which we 
obtain the propagators via the relations above. We employ the diagrammatic Hard Loop 
approach but we have checked that similar expression can be derived from Vlasov transport 
theory [13]. The contribution from the quark loop to the gluon self energy is of the form 



nM^(p) 



tr Yl^S{Q)^''S{K)\ 



(22) 



where S denotes the "bare" quark propagator and Q = K — P. Summing the 11 and 12 
components of the Keldysh representation leads to 



n^(P) = -iNjg' 



d^K 
(2^ 



(go^o + q ■ k) 



Kf{Q)Kr{K) + Ka{Q)Kf{K) 
vKa{Q)Ka{K) + Kr{Q)Kr{K) 



(23) 



where we neglect the fermion mass and write the fermion propagator as Sr^a,f{K) = 
^Ar^a,f{K). The distribution function which appears in the symmetric propagator has 
the form 



k /T 

/^(k) = riFik) - inUkf^i}^ ■ nf + 0(e 



(24) 



with nplp) a Fermi-Dirac function. The last two terms of the integrand vanish after inte- 
gration over ko- Temperature independent terms will be dropped in the following. Shifting 
variables K —K + P in the first term, we found that the first two terms give the same 
contributions to the final result. This is still true for a non-equilibrium distribution which 
satisfies /(k) = /(— k). Then, 



n^(P) = 4vriV^/ y__/^(k) 



{2k^ - Po^ - k ■ p) 



1 



+ (2A;2+pofc-k-p) 



P2 — 2A;po + 2k ■ p — i sgn(/i; — p^e 
1 



P^ + 2fcpo + 2k ■ p — z sgn(— A; — po)^ 



(25) 



Adopting the hard loop approximation, we assume that the internal momenta are of order 
T and therefore much larger than the external momentum which is of order gT [l|. The 
integrand in the square bracket can then be expanded in powers of the coupling and the 
leading term is of the form 



2fc2 



+ 



2fc2 



-2/cpo + 2k ■ p — i e 2A;po + 2k ■ p 



I e 



(26) 
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It can be easily shown that after integrating over dQ, this contribution vanishes. The next 
to leading term comes from the following 4 terms in the expansion of the integrand in the 
square bracket 

-Pok - k ■ p pok-k-p 



-2kpQ + 2k • p — i e 2kpQ + 2k • p + i e 

(27) 



{-2kpo + 2k ■ p - i e)2 {2kpo + 2k ■ p + i e)2 ' 

After integrating over dQ, we find that the first two terms in Eq. (!27|) give the same result; 
the last two terms also contribute equally. As a result, the retarded self energy can be 
expressed as 

^UP) = fkdkdQMk) . (28) 

) (27r)4 J •'^^ ^(k. p + Po±i^)2 ^ ^ 

Here, three-momenta with a hat denote unit vectors. 

Using the distribution function fl2^ it is now straightforward to calculate the retarded 
self energy: 

^RiP) = i-i^f E rk^^i){k)dk f M/(,)(.)ci. , (29) 



with 



^[o){k) = npik), 

pk/Tu 

^(o)(s) 



s + 



sHl-s^) sin2« (1-^2^2 



vl/(i)(s) = cos^ a + ^ ^^^^^ . (30) 



Here, a is the angle between n and p and s = k ■ p. This integral can be performed 
analytically and the leading contribution for the isotropic term is 

n^(o)(P) = f-^Nj rk^^o){k)dk f ^>(^{s)ds = n/-^ f^^^ Po+P + ^e _ \ 
^ ' 2tx^ Jo J-1 6 \2p pn — » + ze / 



To order ~ ^, the result is 

t2 noo p1 



n^(i)(^) = 1^^/ fk<^(imdk J_^^i^{s)ds 

^ ^ (32) 



g'T' (1 cos(2a) , 



2 

cos(2a) -7^(1 + 3 cos(2a)) 



2p2 



In the Hard Loop limit the gluon-loop contributions to the gluon self energy have the same 
structure as the one due to a quark loop. (For the same reason, the gluon self energy in 
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the Hard Loop approximation is gauge invariant.) We can simply replace Nf^-^ by m|, to 
generalize the QED result to QCD. In the following, the Debye mass m|, is to be understood 
as 



This is independent of ^ as all viscous corrections shall be written explicitly. The isotropic 
distribution function f^^o in our case is just a sum of the Fermi and Bose distributions (with 
appropriate prefactors counting the number of degrees of freedom jsl). For Nj massless 
quark flavors and Nc colors, 

^ iNf + 2N,). (34) 



D 
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Analogously, the advanced self energy is given by 



m 



D 



Po_ Po+p-ie _ _^ 
2p po-p-ie . 



(35) 



and 



ni(i)(P) 



m 



D 



'l cos(2a) 
6 ^ 2 



+ ni(o)(P) 



Pi 



cos(2a) - ^ (1 + 3 cos(2a)) 
2p^ 



(36) 



Using Eqs. (1161) and (l20l) . it is straightforward to obtain the temporal component of the 
retarded propagator in Coulomb gauge 



D 



R{0) 



P 



Po+p + te 
m^{ — in ■ 



-1 



2p po- p + ie 



(37) 



+ n^(o)[cos(2a)-|^(l + 3cos(2a))] 



R{1) 



2p^ 



ip' 



m, 



^£0 In 20±P±l£ _ i))2 

< zn nn—n-i-iF ' ' 



(38) 



D\2p po—p+ie 

Similar results can be obtained for the advanced propagator. These results are identical to 
the ones obtained within the transport theory approach. 

Next we calculate 11^ within the hard loop approximation. Eq. f|T4|) shows that this 
quantity is necessary to obtain the resummed symmetric propagator out of equilibrium. 
Summing the 11 and 22 components of the Keldysh representation, we obtain 



n^(P) = -iNfg' 



(2^ 



(go^o + q ■ k) 



Kf{Q)Af{K) - {Kr{Q) - Ka{Q)) 



X {Kr{K) - Ka{K)) 



(39) 



Using Kr{Q) — Ka{Q) = — 27ri sgn(go)5(Q^) and the hard loop approximation, we find that 
the symmetric self energy can be expressed as 



n^(P) = UN,g\' J ^^fMifFik) - I 2 



i2ny 



p 



o{s H ) +d{s ) 

p p . 



(40) 
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Again, we can expand the anisotropic distribution function to order ^ and finally arrive at 



n^(0)(^) = -27rtml'^e{p^-pl) 



^F(i)iP) = ^mml'^ (sm^ a + (3 cos^ a - l)^"! Qip'' ~pI) . (41) 



Next, we can calculate the symmetric propagator. We first consider the isotropic case 
with ^ = and perform a Taylor expansion assuming pq 0: 



which follows from 



hm (In P^±P±11 _ In E^L±Z^) = . (43) 

po^o pq — p + ie Pq — p — ie 

Similarly, when po is small, the distribution function of on-shell thermal gluons is 

2T 

(1 + 2nB) sgn(po) = — • (44) 
Vo 

In the above equation, terms that do not contribute to the symmetric propagator in the 
static limit have been neglected. Finally, we can determine the temporal component of the 
symmetric propagator explicitly in the isotropic limit: 

We now consider the contribution to order C,- From Eq. fl2T|) . the gluon distribution function 
can be expanded as 

fB = /b(o) + = ^ - (46) 

bol 2\po\ 

There are 4 contributions at linear order of ^ as shown in Eq. (1211) . The calculation is similar 
to the isotropic case. In addition, we need the following Taylor expansion to linear order in 
Po- 



m^(l — 3 cos(2a;)) m|, cos(2a; 



+ 



Qp(^p2 _|_ 777,2^)3 2p{p'^ + m'jj^'^ 



Po 



n-,,(P)-n%,(P) = =p(Mp„. (47) 

Now the 0{C,) term of the symmetric propagator in the static limit can be derived from the 
above equations 

n*L / n^ 3mTml . ^ AmTm% 2 

D Fn^(Po = U) = — ^-ttt sm a -— ^— :r(sm a 1 . (48) 
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Heavy quark potential in an anisotropic plasma: In the real time formalism, the static 
heavy quark potential due to one gluon exchange can be determined through the Fourier 
transform of the physical "11" component of the gluon propagator in the static limit: 



V{r,0 = -g'Cp 



(27r)3 



(e*P--l)(Z}*^(po = 0,p,0 



l)^ D 



1 



D* 



D 



(277)3 
rf3p 



-1)-(D*^ + D*^ 

(2.)3(^^^-^-^)^^"- 



(49) 



In the static limit, | [d*^ + D*^ = D*^ = D*\. The Fourier transform of this quantity 
gives the real part of the potential which has been discussed previously in Refs. 13|, ll^l and 
which determines the quarkonium binding energies. Here, we instead consider the imaginary 
part which comes from the Fourier transform of the symmetric propagator. From Eqs. (145|) 
and (HHj) . the isotropic contribution is given by 



Im V(o) (r 



-g'c, 



rf3 



P 



(27r)3 



p {p^ + m\) 



2 \2 



g^CpT 

Ait 



(50) 



with 



2 / dz- 



smizr) 



:z2 + l)2 Zf \ ' ^^^^ 

and f = rmu. This result has been derived before in Refs. js], 0|- The term of order ^ can 
be expressed as 



ImeV^(i)(r) = -g'Cp^ 



d' 



P 



(2vr)^ 



Ap {p^ + m']-)Y 

g^Cp^T 



2 2'7rTmi) ^ 
sm a TTT^ [ sm a — 



p {p^ + m|))3 



47r 



(52) 



where 9 is the angle between r and n and 

3 



1 



9 ^sm Lzr , 9 
sin^e — ' ''"""2, 

zr 



+ (1 -3cos"e)G(f,2) 



(53) 



^2(r,^) 
with 



dz 



^z 

3^ 



(;z2 + 1)3 



1-3 



G(f,;^) 



2 2 /lA sin(2;f) , 

3 J zr 



rzcosyrz) — sm(^rzj 
irzY 



3cos2^)G(f,^) 



(55) 
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FIG. 1: Imaginary part of tlie static potential as a function of distance (r = rrriD)- The vertical 
axis is Im V/{g'^CpT). The curves, from bottom to top, correspond to an anisotropy of ^ = and 
^ = 1, = 0, = 7r/2. 



The result is shown in fig. [TJ The imaginary part decreases with ^ (or with the viscosity, 
respectively). When f is small, we can expand the potential. This is relevant for bound 
states of very heavy quarks whose Bohr radii ~ I/Iq'^Mq) are smaller than the Debye 
length l/rrir)- For the imaginary part, at leading order, the corresponding functions take 
the following forms: 

0(f) = -^f2(-4 + 37£; + 31nf) , 

iP^{r,9) = -i-f2[123 - 907E-901nf + cos(2^)(-31 + 307ij + 301nf)] , 
600 

Mf,0) = ^f2(-4 + 3cos(2e)) , (56) 
where is the Euler-Gamma constant. At leading logarithmic order then 
to V(r,0 = -^f^ 3-cos(2^)\ 

Treating the imaginary part of the potential as a perturbation of the vacuum Coulomb 
potential provides an estimate for the decay width. 



Att J ' ' " f \3 20 J 

-?CVM|l 2j^^^^- ^^^^ 

Here, Mq is the quark mass and \E'(r) is the ground state Coulomb wave function. Thus, at 
leading order in the deviation from equilibrium (i.e., viscosity) the quarkonium decay width 
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is smaller. For a moderate anisotropy ^ ~ 1, F decreases by about 50% as compared to an 
ideal, fully equilibrated plasma. 
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Note added: While this paper was in the final stages of preparation, another paper ap- 
peared where the imaginary part of the heavy-quark potential at non-zero plasma anisotropy 
is being considered |15j. The authors point out that the C(0 correction contributes already 
at leading non-trivial order to the quarkonium decay width, which agrees with our finding. 
Their result differs from ours (numerically) but can be reproduced if the second contribution 
on the r.h.s. of Eq. f[T^ is omitted. From private communication with M. Laine, this appears 
to be due to a different setup of the non-equilibrium system: Ref jl5| considers a situation 
where the soft gluons are in equilibrium at a temperature T while the hard gluons are out of 
equilibrium and are characterized by a different hard scale T'. In contrast, we assume that 
the deviation from equilibrium follows viscous hydrodynamics, Eq. ([2]). 
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